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ABSTRACT 

Lyman-a absorption is a prominent cosmological tool for probing both galactic halos and the intergalactic 
medium at high redshift. We consider a variety of sources that can be used as the Lyman-a emitters for this 
purpose. Among these sources, we argue that quasars are the best probes of the evolution of massive halos, 
while 7-ray bursts represent the cleanest sources for studying the reionization of the intergalactic medium. 
Subject headings: galaxies: high-redshift, cosmology: theory, galaxies: formation 



1. INTRODUCTION 

Over the past three decades, quasars have been tradition- 
ally used to probe the intervening intergalactic medium (IGM) 
at redshifts z < 6 when this medium was highly ionized. 
The polarization anisotropies of the microwave background 
measured by WMAP indicate that the IGM was signif icantly 
ionize d at a much higher redshift z ~ 17 ± 4 ( Spergel et al. 
2003), possibly implying two separate reionization phases 
separat ed by a p rolonged period of nearly c omplete recombi- 
nation dWvithe & Loebl2003alb1:ICenl2003T) . Attention is now 
shifting towards direct observations of the epoch of reion- 
ization. The process of reionization started with the emer- 
gence of individual H II regions centered on discr ete sources, 
which e ventually overlapped [se e reviews bv lBarkana & Loea 
(2001); Loeb & Bar kand feOOlM . Observationally, a central 
goal is to determine the average neutral fraction of the IGM, 
and possibly even map the topology of the H II islands, as a 
function of redshift. The main problem with using resonant 
Lya absorption by neutral hydrogen for this purpose is that 
the corresponding resonant cross-section is too large; a neu- 
tral fraction as small as ~ 10~ 3 is already sufficient to suppress 
the quasar flux below detectability. For example, the spec- 
trum of the Sloan Digital Sky Survey (SPS S) quasar at z = 
6.28 sh ows complete absorption [the so-called |Gunn-Petersonl 
i ll 9651) trough] of its spectral flux shortw ard of the Lya reso - 
nance wavelength at the quasar redshift ( Bec ker et alJ UOQl). 
but it remains unclear whether the IGM is highly ionized 
or whether it contains poc kets of neutral gas at this redshift 
(Barkana 2002: lFan et all2002h . 

The use of other atomic lines is problematic for o ther 
reaso ns. For example, absorption by metal lines dOhl 
l200l can only probe regions that had been previously 
metal-enriched by supernovae and the ionization state of 
this disturbed gas may not be representative of the entire 
IGM (since most metals may reside within H II bubbles 
near their source galaxies). Another possible probe, 21 
cm absorption by neutral hydrogen, suffers from a cross- 
section that is too small (the opposite problem from Lya) 
leading to a very low optical depth and requiring the use 
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of very bright radio sources (> 10 mJy) that should be 
rare in the early univer se dCarilli. Gnedin. & Owed l2002t 
Furla netto & Loebl 120021) . Detection of 21 cm absorption 
or emission against the microwave b ackgrou nd is tech- 
nolog ically challenging Q bzzi, Madau, Meiksin, & Rees 
2000; llliev. Shapiro. Ferrara. & Martell l2002t 

Furla netto. Sokasian^c^Hernquistl 120031) and 

is likely to be overwhelmed by fo regrounds 
dPi Matteo. Perna. Abel. & Reesl200llLoebll99'6l) . 

It is therefore important to explore any possible way of 
using the Lya absorption signature near the source in or- 
der to measure the IGM neutral fraction. In fact, for pho- 
ton wavelengths longer than the Lya resonance at the source 
redshift, the absorption cross-section is weaker (since these 
photons are offset from resonance) and so the quasar spec- 
trum should gradually recover its unabsorbed spectral flux 
level as the wavelength becomes redder. The detailed spec- 
tral shape of the red damping wing of the Gunn-Peterson 
trough can potentially be used to infer the neutral hydro- 
gen fraction in the I GM during the epoch of reionization. 
Miralda-Escude ( 1998) calculated this shape for the idealized 
case of a weak point source embedded in a neutral uniformly 
expanding IGM. In reality, however, the ionizing effect of a 
quasa r on the surrou nding IGM [the so-called "proximity ef- 
fect" dBaitlik. Puncan. & Ostrik eJ 19881)1 complicates thi s in- 
terpretation dCen & Haimanl l2000: Madau & Reesll2000l) . In 
addition, cosmological infall changes the IGM density and ve- 
locity field near bright early quasars because these ar e likely 
to be hosted by a massive galaxy; Barkana & Loeb (2003) 
have shown that this effect can be used as an important tool in 
itself, as a way to measure directly the masses of the host dark 
matter halos of quasars, and to test fundamental aspects of the 
processes of galaxy and quasar formation. Indeed, we explore 
this tool further in this paper and consider how infall around 
quasars can be probed over a broad range of halo masses and 
redshifts. 

Given that quasars and their massive hosts perturb the IGM 
around them both radiatively and gravitationally, the follow- 
ing question arises: Are there better sources than quasars 
for probing the Lya damping wing during reionization? Our 
answer is positive. These sources are the afterglows of 
Gamma-Ray Bursts (GRBs), suggested before for this pur- 
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pose ( Mira lda-EscudelH998t lLoehll200l but with additional 
advantages that we highlight in this paper. 

GRBs are the brightest electromagnetic explosions in the 
universe, and should be detectabl e out to r edshifts z > 10 
(Lamb & Reicharl2000t ICiardi & Loebl2000l) . High-redshift 
GRBs can be easily identified through infrared photometry, 
based on the Lya break induced by absorption of their spec- 
trum at wavelengths below 1.216/.tm[(l +z)/10]. Follow- 
up spectroscopy of high-redshift candidates can then be per- 
formed on a 10-meter-class telescope. There are three main 
advantages of GRBs relative to quasars: 

• The afterglow flux at a given observed time lag after 
the 7-ray trigger is not expected to fade significantly 
with increasing redshift, since higher redshifts trans- 
late to earlier times in the source fram e, during which 
the af t erglow is intrinsically brighter (Lamb & Reichart 
2000; Ciardi & Toebl 120001) . For standard afterglow 
lightcurves and spectra, the increase in the luminosity 
distance with redshift is compensated by this cosmic 
time-stretching effect. 

• In the standard ACDM cosmology, galaxies form hi- 
erarchically, starting from small masses and increasing 
their average mass with cosmic time. Hence, the char- 
acteristic mass of quasar black holes and the total stellar 
mass of a galaxy were smaller at higher redshifts, mak- 
ing th ese sources intrinsically fainter ( Wvithe & Loeb 
2002). However, GRBs are believed to originate from 
a stellar mass progenitor and so the intrinsic luminos- 
ity of their engine should not depend on the mass of 
their host galaxy. At high redshifts, a GRB should 
strongly dominate the luminosity of its host galaxy. 
Hence GRBs are expected to be increasingly brighter 
than any other competing source as higher redshifts are 
considered. 

• Since the GRB progenitors are believed to be stellar, 
they likely originate in the most common star-forming 
galaxies at a given redshift rather than in the most mas- 
sive host galaxies, as is the case for bright quasars. Low 
mass host galaxies induce only a weak ionization ef- 
fect on the surrounding IGM and do not greatly perturb 
the Hubble flow around them. Hence, the Lya damp- 
ing wing shoul d be closer to the idealized unperturbed 
IGM case ( Miralda-Escude 1998) and its detailed spec- 
tral shape should be easier to interpret. Note also that 
unlike the case of a quasar, a GRB afterglow can itself 
ionize at most ~4x \0 a Es\Mq of hydrogen if its UV 
energy is E51 in units of 10 51 ergs, and so it should have 
a negligible cosmic effect on the surrounding IGM. 

Although the nature of the central engine that powers 
the relativistic jets of GRBs is unknown, recent evidence 
indica tes that GRBs trace the formation o f massive stars 
JBloom et alJ l20()l iKulkarni et alJ pOOCt iTotanil fl997t 
IWiiers et alJ 119981 iBlain & Nataraianl 120001) . Since the 
first stars are predicte d to b e pre dominant ly massive 
iAbel. Brvan. & Normanl 120021 iBromm. Coppi. & Larson! 
120021) . their death might give rise to large numbers of 
GRBs at high redshifts. Detection of high-redshift GRBs 
could probe the earliest epochs of star formation, one 
massive star at a time. The upcoming Swift satellite (see 
http://swift.gsfc.nasa.gov/), scheduled for launch in May of 
2004, is expected to detect about a hundred GRBs per year. 



IBromm & Loebl J2002I) calculated the expected redshift dis- 
tribution of GRBs; under the assumption that the GRB rate is 
simply proportional to the star formation rate they found that 
about a quarter of all GR Bs detected by Swif t should originate 
at z > 5 (see also Lamb & Reichart 2000). This estimate 
is rather uncertain because of the poorly determined GRB 
luminosity function. We caution further that in principle, the 
rate of high-redshift GRBs may be significantly suppressed 
if the early massive stars fail to launch a relativistic outflow. 
This is possible, since metal-free stars may experience 
negligible mass loss before exploding as a supernova. They 
would then retain their massive hydrogen envelope, and any 
relativistic jet mig ht be quenched before escaping the star 
jHeger et alJ I2003I) . However, localized metal enrichment 
is expected to occur rapidly (on a timescale much shorter 
than the age of the then-young universe) due to starbursts 
in the first galaxies and so even the second generation of 
star formation could occur in an interstellar medium with 
a significant metal content, resulting in massive stars that 
resemble more closely the counterparts of low-redshift GRB 
progenitors. 

In this paper we compare quantitatively the spectral profile 
of Lya absorption for quasars and GRB afterglows, including 
both resonant absorption and the damping wing and taking 
into account the ionization and infall effects induced by their 
host galaxies and halos on the surrounding IGM. We present 
our detailed model in §2, and then illustrate it in §3. In §4 
we show how quasars can be used as a tool for measuring 
infall, which probes the properties of massive halos and their 
formation process. In §5 we show how, on the other hand, 
GRBs can be used to cleanly probe the reionization of the 
IGM. Finally, we give our conclusions in §6. 

2. MODELING Lya TRANSMISSION 
2. 1 . Summary of Simple Models 

We assume throughout a standard ACDM cosmology. For 
the contributions to the energy density, we assume ratios rel- 
ative to the critical density of fl m = 0.27, Oa = 0.73, and 
fib = 0.044, for matter, vacuum (cosmological constant), and 
baryons, respectively. We also assume a Hubble constant 
Hq = 7 1 km s -1 Mpc -1 , and a primordial nearly scale invari- 
ant power spectrum (with a slowly varying spectral index) 
with erg = 0.84, where erg is the root-mean-square amplitude 
of mass fluctuations in spheres of radius 8 h~ l Mpc. These 
parameter values are based on recent temperature anisotropy 
measurements of the cosmic microwave background in com- 
bination with other large scale stru cture measurements, and 
the full set of parameters is given by Sperg el et alJ (|2P03). 

An ionizing source embedded within the neutral IGM can 
ionize a region of maximum proper size 
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assuming that recombinations are negligible and that all A^t 
ionizing photons are absorbed by hydrogen atoms. More 
generally, the evolution of an H II region, including a non- 
steady ionizing source, recombinations, and cosmological ex- 
pansion, is given by the e xpanding Stromgren sphere equation 
( Shapi ro & Girouxl 19871) . 
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where V is the comoving volume of the H II region, dN 1 / dt 
is the number of photons per unit time output by the source, 
n° H is the z = cosmic mean number density of hydrogen, and 
the case B recombination coefficient for hydrogen at T = 10 4 
Kisa B = 2.6x 10" 13 cm 3 s" 1 

When dealing with Lya emission and absorption, it is use- 
ful to associate a redshift z with a photon of observed wave- 
length A b s according to A D b s = A Q (1 +z), where the Lya wave- 
length is X a = 1215. 67A. In the absence of infall, the optical 
depth Tdamp due to the IGM damping wing of neutral gas be- 
tween Z] and z? (where z\ < z ? < z) can be calculated analyt- 
ically ( Miralda-Escude 1998). For a source above the reion- 
ization redshift, zi is the redshift corresponding to the blue 
edge of the H II region, and zi is the reionization redshift. In 
the limit where zi is very close to z while z\ is not, the full 
formula can be approximated as 
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where AA = A a (z—Zz) and in this limit the result is indepen- 
dent of Zl • 

In addition to the damping wing due to neutral gas that lies 
outside the H II region, photons may also be absorbed near the 
Lya resonance by the small neutral hydrogen fraction within 
the ionized region. In the absence of clumping, this resonant 
absorption produces an optical depth 



Tres = TCpXm A 



(4) 



assuming gas at relative densit y A = p„/p„ and a neut ral frac- 
tion xri, where the standard dGunn-Petersonl H%5l) optical 
depth is 
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The neutral fraction is given by ionization equilibrium. The 
optical depth at a proper distance R from a source that spews 
hydrogen-ionizing photons into the IGM at a rate dN 1 /dt is 



r res =0.731 A 2 
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where 5# is the frequency-averaged photoionization cross 
section of hydrogen (which depends on the ionizing photon 
spectrum), and this result uses the value of a# given above. 

Previous stud ies that included the absorption of gas in the 
H n region (ICen & Haimanl2000t iMadau & Reesll2000l) have 
all assumed that the mean gas density in this region is equal 
to the cosmic mean density, and for the most part have treated 
gas clumping very simply, taking a constant clumping fac- 
tor C =< A 2 > / < A > 2 — 10 to enhance both r res and the 
recombination rate du ring the growth of the H II region [see 
eq- ©] ■ More recently, Haiman & Cen ( 2002) used numerical 
simulations directly, while lHahnanT(l2"002l) assumed a lognor- 
mal distribution of gas densities. In the following section we 
account for infall and also use a realistic redshift-dependent 
distribution of gas clumping that is based on numerical simu- 
lations. 



2.2. Gas Infall and Clumping 

In the standard CDM-dominated picture of galaxy forma- 
tion, galaxies form within dark matter halos that themselves 
form only in regions of high initial density. Since the den- 
sity field is generated according to an initial power spectrum, 
the region around a halo is strongly correlated with the region 
that leads to the halo itself and the former also tends to have 
a higher than average initial density. As the dense central re- 
gion undergoes gravitational collapse, the surrounding region 
is drawn in towards the forming halo. Thus, by the time a 
given dark matter halo has finally virialized, a much larger 
region of surrounding gas has already acquired a significant 
infall velocity and, by mass conservation, a large overdensity. 
Infall affects resonant Lya absorption through A [eq. (0}] and 
also affects the damping wing, for which the optical depth is 
no longer given by a simple formula such as eq. Q but must 
b e numerically in tegrated over the infall profile. 

Barkana (2004) developed a model for calculating the ini- 
tial density profile around overdensities that later collapse to 
form virialized halos. This model, which we adopt in this pa- 
per, accounts for the so-called "cloud-in-cloud" problem: if 
the average initial density on a certain scale is high enough 
to collapse by a given redshift, this region may be contained 
within a larger region that is itself also dense enough to col- 
lapse; in this case the original region should be counted as 
belonging to the halo with mass corresp onding to the larger 
collapsed region. Technically, Barkana] (120041) adopted the 
definition of hal os in t erms o f the initial density field as pro- 
posed by Bond et all (119911) . who improved on the original 
Press & Schechteg (119741) model by including the cloud-in- 
cl oud prob l em in the definition and in the statistics of halos. 

Barkana (2004) found that the mean initial density profile 
around an overdensity that leads to a halo depends on both the 
mass of the halo and its formation redshift. At high redshift, 
low mass halos tend to have a stronger density enhancement 
around them that extends out to larger distances (in units of 
the initial radius containing the virial mass of the eventual 
halo); this is a result of the presence of stronger correlations 
on smaller scales in the standard cosmological model. How- 
ever, at lower redshifts, low mass halos tend to be surrounded 
by initial regions of relatively low density. This can be un- 
derstood as follows: low mass halos that were initially sur- 
rounded by a high overdensity would most likely end up as 
part of a larger halo, since even large halos are common at 
low redshift; therefore, isolated low mass halos at low red- 
shift should typically be surrounded by initial regions of low 
density, preventing them from becoming part of a larger col- 
lapsed halo. 

Starting from this mean initial profile, spherical collapse 
can be used to obtain the final density and velocity profile 
surrounding the virialized halo. In reality, there will be some 
variation in the initial profiles, with each lead ing to a d ifferent 
final profile. This scatter was also studied by Bar kana! J2004I) . 
but in this paper we adopt the mean profile described above, 
which accounts for the overall trends in the infall profile as a 
function of the halo parameters. 

Once the dark matter infall profile is fixed, we assume that 
the infalling gas follows the dark matter at large distances 
from the halo, where pressure gradients are negligible com- 
pared to gravitational forces (i.e., we only consider halos well 
above the Jeans mass of the intergalactic gas). When the 
gas finally does fall into the halo, incoming streams from 
all directions strike each other at supersonic speeds, creat- 
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ing a strong shock wave. Three-dimensional hydrodynamic 
simulations show that the most massive halos at any time in 
the universe are i ndeed surrounded by strong, quasi-spherical 
accretion shocks (Abel, Brvan, & Norman 2002; Keshet et al. 
2003). These simulations show that along different lines of 
sight, the shock radius lies at a distance from the halo center 
of around 1-1.3 times the halo virial radius. We adopt 1.15 
as an illustrative value, noting that the results are not altered 
substantially as long as the shock radius is fairly close to the 
virial radius. 

We calculate gas infall down to the radius of the accretion 
shock, and neglect any Lya absorption due to the post-shock 
gas. The post-shock gas is heated to around the virial temper- 
ature which, in halos we consider, is > 10 4 K. In halos where 
the temperature is far greater than 10 4 K the gas is highly 
collisionally-ionized in addition to the effect of the intense 
radiation field near the source, and the absorption should be 
weak. If the gas is instead shocked to a temperature close to 
10 4 K then it will subsequently cool and likely collapse onto 
the galactic disk. Hydrodynamic simulations are required in 
order to check for the possibility of having a thin cold shell 
of shocked gas that may produce some absorption, but in any 
case such absorption will be seen well to the blue side of the 
absorption due to the infalling pre-shock gas. 

In addition to the increase in the mean gas density due to 
infall, the gas is also clumped, and this has an additional ef- 
fect on the recombination rate and on Lya absorption. For 
the distribution of gas clumps we adopt the distribution that 
Miralda-Escude et al. (2000) constructed based on numerical 
simulations at z = 2-A and extrapolated to other redshifts: 
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where clumps at an overdensity Ac acco unt for a volume frac- 
tion Py (Ac)dAc of the gas. Following lMiralda-Escude et alJ 
(2000), we adopt (3 = 2.5 for z > 6 and slightly lower values 
of j3 at lower redshifts, set 6q = 7.61/(1 +z), and fix the other 
parameters by requiring t he total volume and mass fr actions 
to be normalized to unity. Miralda-Escude et al. ( 20001 com- 
pared their clumping distribution only to post-reionization 
simulations, and only up to Ac ~ 100; however, we are justi- 
fied in adopting this distribution (up to A max = 100) since we 
apply it only within H II regions, where the gas has already 
been photoionized and heated, and we also demonstrate be- 
low that all our observable predictions are insensitive to the 
high-density tail of the distribution. In our model we assume 
that eq. Q applies with Ac interpreted as the extra overden- 
sity due to clumping, measured relative to the mean gas den- 
sity A (which itself may be higher than the cosmic mean, due 

to infall). 

Miralda-Escude et al. (2000) modeled reionization with a 
model that assumes that all clumps up to some critical Ac 
are fully ionized while clumps above this Ac are completely 
neutral. This model is inadequate for our purposes, since the 
small neutral fraction expected in the low-density clumps is in 
fact crucial in determining the total Lya absorption, and in ad- 
dition high-density clumps are in reality not self-shielded and 
fully neutral except at very high densities that do not affect 
our results. We thus assume that clumps are optically thin in- 
stead, and find the neutral fractionxei( Ac) separately for each 
Ac based on ionization equilibrium of gas at that overdensity. 
This implies that the neutral fraction essentially increases in 
proportion to Ac (except when ;tHi(Ac) — > 1, which occurs 



only for extremely high overdensities and does not affect our 
results). Thus, we may define one type of effective clumping 
factor in terms of the total neutral fraction compared to the 
value it would have in the absence of clumping: 
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where A is the enhancement of the mean density due to gas 
infall. 

The recombination rate of a given gas element and the op- 
tical depth to ionizing photons are both proportional to C Xm . 
However, Lya absorption is not determined simply by the to- 
tal number density of hydrogen atoms, since it results from 
resonant absorption by separate gas elements at a variety of 
overdensities; even if high-density clumps produce complete 
absorption, photons at wavelengths that resonate with gas in 
voids are strongly transmitted. Cosmological volume is dom- 
inated by voids, and furthermore the mean fraction of a given 
line of sight covered by gas of a given overdensity is simply 
equal to the volume fraction occupied by gas at that overden- 
sity, regardless of the topology of the c lumps [e.g., see foot- 
note in §2.5 of Barkana & Loeb (2002)]. Thus, the effective 
clumping factor for Lya transmission C]_ ya is determined by 



-[r G pC LyQ Aini(A)] _ 
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P v (A c )dA c . 

(9) 

Note that when we use C]_ ya , we are in effect averaging over 
the density fluctuations in the IGM. This could in principle be 
compared directly to observations averaged over many differ- 
ent lines of sight through the IGM. Any particular line of sight 
is expected to show significant fluctuations (in the resonant 
absorption component) relative to the average pattern that we 
predict. These fluctuations are the high-redshift equivalent 
of the Lya forest of absorption features seen in the IGM at 
lower redshifts. Note that unlike resonant absorption, damp- 
ing wing absorption is determined by the cumulative effect of 
a large column of gas (typically of order 1 Mpc long) and is 
therefore unaffected by fluctuations and clumping that occur 
on much smaller spatial scales. 

The mean free path to ionizing photons depends on the total 
neutral fraction and thus on C Xm , since absorption of ionizing 
photons is not a resonant process but is instead cumulative 
along a line of sight. While C]_ ya is typically dominated by 
low values of Ac < 1, C Xm is sensitive to the high integration 
limit A max . However, this sensitivity is a mathematical arti- 
fact, not a real physical dependence of the mean free path of 
ionizing photons; it shows that if gas clumps were uniformly 
distributed, so that the one-point distribution of gas density 
in arbitrarily small volumes were given by eq. Q, then the 
high-density tail of clumps would indeed dominate the value 
of C Xm . However, in reality gas clumps in the ionized IGM at 
densities > 20 are found only around or inside localized, rare, 
virialized halos. The clumping distribution is only correct as 
an average over large volumes that contain a large statistical 
sample of voids, filaments, and halos [ but no te that the simu- 
lations fitted by Miralda-Escude et al. (2000) could not probe 
a possible population of minihalos]. An ionizing photon, for 
example, may travel a long distance through a void before 
encountering any clump of significant overdensity, and thus 
the mean free path of ionizing photons is determined by the 
topology of the clu mps and not just by th e overa ll one-point 
distribution of Ac- Miralda-Escu de et alJ (|2000) solved this 
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problem approximately, finding a rough formula for the typ- 
ical distance d (defined as a mean free path) along a random 
line of sight between encounters of clumps with Ac > Ac: 



t/(A c )«60kms- 1 i/- 1 (z) 
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where H(z) is the Hubble constant at redshift z. If, for a given 
Ac, we consider a line of sight of length / through an H II 
region, then the probability of encountering any gas clump 
with A c > A c somewhere along this line of sight is 
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In this paper, we consider absorption along relatively short 
lines of sight from a source. Most such lines of sight will 
not contain any clumps with very high density A c . Thus, we 
fix our maximum A max so that a line of sight of the length 
we are considering is reasonably likely to go through at least 
one clu mp of that high a density. Since Mira lda-Escude et alJ 
(2000) compared eq. J10I to simulations averaged over the 
mean IGM, in regions overdense due to infall we replace the 
line of sight length / by an effective length weighted by the 
overdensity A. In our model we set A max according to P any (> 
A max ) = 0.5 (i.e., a typical line of sight), but the absorption 
profiles are only weakly sensitive to the adopted probability. 
The resulting values of A max range from ~1 for the smallest 
(~ 0.1 Mpc) H II regions up to ~ 50 for the largest (~ 10 
Mpc) H II regions. 

2.3. Basic Halo Parameters 

The absorption profile due to H I in the IGM depends on 
several basic parameters. The first is the source redshift zs- 
The second is the source halo mass M, which determines the 
density and velocity profile of infalling gas, as well as the 
virial radius and thus the shock radius. The remaining three 
describe the source ionizing photon production. One parame- 
ter is dNry/dt, the total rate at which hydrogen ionizing pho- 
tons from the source enter the IGM. The second parameter is 
the age of the source ts, which is the period of time during 
which the source has been active (with an assumed constant 
dNry/dt, for simplicity). Note that the source lifetime is only 
important prior to the end of reionization, when the lifetime 
affects the size of the surrounding H II region. The final pa- 
rameter is an, the frequency-averaged photoionization cross 
section of hydrogen, which affects the absorption and depends 
on the spectrum of the ionizing photons. 

Bright quasars are powerful ionizing sources and so the 
stellar emission of their galactic host can be neglected. In or- 
der to predict the ionizing intensity for a quasar with a given 
continuum flux, we assume a typical quasar continuum spec- 
trum as follows. We adopt a power-law shape of F v oc v~~° M 
in the rest-frame range 1 190-5000A based on the SDSS com- 
posite spectrum (IVanden Berk et al.l2001l) . and F v oc v~ 1S1 at 
500-1 190A using the co mposite quasar spe ctrum from the 
Hubble Space Telescope ( Telfej^eUyJ|2002|) . Based on ob- 
servations in soft X-rays ( Yu an et aLlll9 98). we extend this 
power-law towards short wavelengths. We assume that the 
brightest quasars shine at their Eddington luminosity (but we 
also compare to the case of a tenth of Eddi ngton), and we note 
that for the SDSS composite spectrum ( Vande nBerk et ail 
1200 ll) . the total luminosity above 1190A equals 1.6 times 



the total continuum luminosity at 1190-5000A. Thus, ion- 
izing photons stream out of the quasar's galactic host at the 
rate N = 1.04 x 10 56 M 8 s" 1 , where M 8 is the black hole mass 
Mbh in units of 10 8 M Q . With these assumptions, the black 
hole mass can be inferred from the observed continuum at 
1350Awiththe conversion F v = 1.74 x 10 30 M 8 erg s" 1 Hz" 1 . 

In order to predict the Lya absorption around a quasar we 
must estimate the mass of its host halo. A tight correla- 
tion has been measured in local galaxies between the mass 
of the central black hole and the bulge velocity dispersion 
jFerrarese & Merrittll2000l iTremaine et afll2002l) . This rela- 
tion also fits all existing data on the luminosity func tion of 
high-redshift quasars within a simple model ( Wvithe & Loeb 
120021 12003d) in which quasar emission is assumed to be 
triggered by merg ers during h ierarchic al galaxy formation. 
We use the best-fit (Femygs e & Merriti2000HTremaine et alJ 
l2002HWvifhe & Loebl2003cl) relation, in which the black hole 
mass Mbh is related to the circular velocity V c at the halo virial 
radius by 



M BH = 8.8 x 10 7 



300 km s" 



M 



(12) 



In general, we assume a source age tg ~ 10 7 — 10 8 yr [see 
(Wvithe & Loeb 2003c) for a discussion]. 

When we consider GRB afterglows at redshifts z > 7, which 
are typically found in relatively small galaxies, the above scal- 
ings imply that the contribution of a central black hole in 
these galaxies is negligible compared to the stars, and we must 
therefore include the stellar output of ionizing radiation. For 
stars, we find it useful to express dN^/dt, for given values of 
M and tg, as follows: 



dN 7 _ m a h „ 10Ji 



dt m p n 



M flh Nja 

ts 



(13) 



where m p is the proton mass, and Ni on gives the overall num- 
ber of ionizing photons per baryon. To derive the expected 
range of possible values of Mon> we note that it is determined 
by 



N ion ~~ ^7 f* Jese 



(14) 



where we assume that baryons are incorporated into stars with 
an efficiency of /*, that iV 7 ionizing photons are produced 
per baryon in stars, and that the escape fraction (out to the 
virialization shock radius) for the resulting ionizing radiation 
is /esc- We adopt /* = 10% (based on a rough comparison 
of models to the observed cosmic star formation rate at low 
redshift), and consider / esc in the range 10%-90%. 

The remaining variable N 7 depends on the stellar initial 
mass function (IMF). We consider two examples of possi- 
ble IMFs.TTie first is the locally-measured "normal" IMF 
of lScaiol (1998), along with a metal licity equal to 1 /20 o f 
the solar value, which yields [using Leitherer et al. (1999)] 
N 7 = 4300. The second is an extreme Pop III IMF, assumed to 
consist entirely of zero metallicity M > 100M Q stars, which 
yields N 7 = 44000 based on the ionization rate per star, the 
stellar spectrum, and the main-sequ ence lifetime of these stars 
llBromm. Kudritzki. & Loebl l2001). In general, though, the 
ionization state of hydrogen depends only on the ionizing 
photons that are actually absorbed by hydrogen; therefore, 
the effective iV 7 for hydrogen, as well as the effective cross- 
section oh [see eq. (|6Jl], are both affected by the presence of 
helium, as explained in the following section. 
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2.4. The State of Helium During Hydrogen Reionization 

Along with the H II ionization front, we calculate the evo- 
lution of the He II and He III fronts with equations analogous 
to eq. ©. Initially, when the IGM is neutral, hydrogen atoms 
absorb all 13.6-24.59 eV photons and helium atoms absorb all 
E > 24.59 eV photons. Since the number density of helium is 
only fue = 7.9% that of hydrogen (assuming 24% helium by 
mass), the He II front easily keeps up with the expanding H II 
front, since most realistic sources produce E > 24.59 eV pho- 
tons relative to 13.6-24.59 eV photons in a much larger ratio 
than fue- However, the He III front may either lag behind (for 
sources which emit relatively few of the E > 54.42 eV pho- 
tons that are needed to fully ionize helium) or may catch up 
with the other two fronts. 

Our calculations of the sizes of the various regions assume 
the same infall overdensity and clumping factor for helium 
and for hydrogen, and numerically solve for the expanding 
Stromgren spheres. However, for most sources the effect of 
helium can be accurately described by one of two limiting 
cases. In the first case (applicable to normal stars), the He III 
front lags far behind and covers a negligible volume compared 
to the common H II - He II front. In this case, all helium 
atoms within the H II region are ionized exactly once, and 
therefore the effective A^ ad for determining the radius of the 
H II region is equal to the total number of ionizing photons 
times 1/(1 +/hc)- In this case as well, all E > 54.42 eV pho- 
tons are absorbed by He II when they reach the He III front, 
so the effective an for hydrogen as well as the effective N™ 
for determining the ionizing intensity are both determined by 
the spectrum and number of photons at E < 54.42 eV; how- 
ever, since this case is by definition accurate only for sources 
which emit very few photons above 54.42 eV, the spectrum 
and intensity are essentially equal to the result obtained from 
all E > 13.6 eV photons emitted by the source. 

In the second case (applicable to Pop III stars as well as to 
quasars), the He III front catches up with the He II front and 
creates a common H II - He III front beyond which both he- 
lium and hydrogen are fully neutral. In this case, all helium 
atoms within the H II region are ionized exactly twice, and 
therefore the effective N™ A is equal to the total number of ion- 
izing photons times 1 /(l +2/# e ). In this case, all E > 13.6 eV 
photons emitted by the source reach the gas within the H II 
region and determine both cr# and A^ on . Although our numer- 
ical results below account for the full complexity, where we 
now give representative numbers for various sources we as- 
sume one of the two limiting cases above, as appropriate for 
each type of source. 

Thus, we find the following cross-sections for the sources 
discussed in the previous subsection: on = 2.8 x 10~ 18 cm 2 for 
stars with a normal IMF, <r/; = 1.2x 10~ 18 cm 2 for stars with a 
Pop III IMF, and a H = 2.3 x 10~ 18 cm 2 for quasars. Note that 
we have neglected some other radiative transfer effects; for 
example, the degraded photons that are re-emitted by helium 
may ionize hydrogen atoms as well, but only a minor fraction 
of the helium atoms are expected to recombine on the short 
source timescale ft. 

The source ionizing spectrum and the presence of helium 
also affect the gas temperature, which depends on photo- 
heating as well as adiabatic cooling, Comp ton cooling, and 
various atomic cooling processes. Abel & Haehnelt ( 1999) 
showed that, beyond heating and cooling processes, the tem- 
perature profile inside an H II region depends on subtle ra- 
diative transfer effects; guided by their results, we simply 



adopt characteristic temperatures of 40,000 K within freshly- 
reionized He III regions and 15,000 K in regions where he- 
lium is not fully ionized. These temperatures affect the case B 
recombination coe fficients of the various ato mic species [we 
adopt values from Verner & Ferland (1996)], which in turn 
affect both the growth of the various ionized regions and the 
final T res at each distance. The final effect of helium that we 
include as well is the enhancement of the recombination rate 
due to the extra electrons from ionized helium, where the en- 
hancement factor equals 1 + fn e in He II regions and 1 + 2/# e 
in He III regions. 

2.5. Transmission After the End of Reionization 

Once reionization of the universe is complete, the mean 
free path of ionizing photons begins to rise. A UV back- 
ground is established which gradually becomes more uniform 
as each point in the IGM sees an increasing number of ion- 
izing sources. In this post-reionization era, we assume that 
the mean free path of ionizing photons is much longer than 
the line of sight that we consider, and thus it is unne cessary 
to estimate the mean free path as was done in 32.21 After 
reionization, there is also no limit to the H II region surround- 
ing a source. In this case, there would still be a limit to the 
He III region around a quasar, if the cosmic UV background 
is not sufficiently hard to fully ionize helium throughout the 
universe. We neglect this, however, since the observable ef- 
fect on Lya absorption would be minor, given that calcula- 
tions with radiative transfer show that the double ionization 
of helium produces a characteristic gas temperature of only 
~ 1 .5 x 10 4 K in regions th at had already been reio nized by a 
softer ionizing background (Abel & Haehnelt 1999). 

In our post-reionization calculations, we add the mean UV 
background to the ionizing intensity of the central source. 
This background determines the average transmission level 
at large distances from the particular galaxy being consid- 
ered. We set the background level at different redshifts so 
that our assumed distribution of gas clumping matches obser- 
vations; specifically, we use eq. (|9} with A = 1 and A max = 100 
and compare it to the me an transmission in the Lyq forest at 
various redshif ts; we use iBernardi et alJ J2003T) for z = 2.5- 
4 and|Songaila & Cowie (2002) for z = 4.5-6, we assume no 
UV background at z > 6.5, and we smoothly interpolate in 
the redshift intervals 4^1.5 and 6-6.5. Note that if much 
higher redshifts (z > 15) are considered, the contribution of 
a cosmic UV background might have to be considered around 
the time of a previous reionization phase as suggested by the 
large-a ngle polarization of the cosmic microwave background 
(Spergel et al. 2003). However, in the examples in this paper 
we refer only to the final phase of reionization at z ~ 7, as- 
suming that the universe had managed to recombine nearly 
completely after a possible earlier reionization phase (see § 1). 

2.6. Absorption by the Host Galaxy 

The galactic host of any potential background light source 
may itself absorb some flux near Lya and thus muddle the 
interpretation of any detected absorption, which is no longer 
exclusively due to H I in the IGM. The two main effects are 
absorption due to H I gas in the galactic host and extinction 
by dust. 

As noted above, absorption due to the damping wing of 
H I in the IGM produces an optical depth that varies approx- 
imately inversely with AA. However, since H I gas in the 
galaxy is concentrated in a compact region compared to the 
apparent distance corresponding to the relevant range of AA, 
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it produces for photons with A D b s = X a ( 1 + z) = X a ( 1 + zs) + A A 
an optical depth 

r — ^(lo^X 1 ?) (^) (Ix) ' 

(15) 

where Nm is the total H I column density in the host galaxy at 
redshift zs- Note that this formula does not assume \z— Zs\ <C 
1 . The inverse-square dependence of Td amp due to the galactic 
host may in principle allow the isolation of this effect from 
damped IGM absorption. However, a column density > 10 21 
cm would obscure much of the region where damped absorp- 
tion could otherwise be measured, and could be partly degen- 
erate with IGM absorption if only a narrow wavelength win- 
dow is available for separating the two effects. As shown by 
iMirald a-Escudei ( I1998I) . such a partial degeneracy may ham- 
per an accurate measurement of the IGM density through the 
normalization of Td amp due to the IGM; however, as we show 
in |5] this degeneracy does not preclude a definite detection 
of absorption due to a neutral IGM, and it mostly disappears 
if the source redshift can be measured accurately from addi- 
tional absorption lines (other than Lya) associated with the 
host galaxy. Note also that in cases where IGM absorption is 
not complete in the Lya blue wing, a neutral gas column in 
the host galaxy could produce an observable, symmetric blue 
damping wing. 

The second source of nuisance absorption is extinction by 
dust in the host galaxy. As a guide to the expected extinc- 
tion we adopt the mean G alactic extinction curve as mod- 
eled by Fitznat nckl ( 119991) . and normalize according to the 
dust-to-gas ratio of the di ffuse Galactic interstellar medium 
iBohlin. Savage. & Drakdfl978l) . We find that extinction at 
Lya is higher than in the V band according to A(X a ) = 
3.3A(V), where A(A) is the extinction (in magnitudes) at 
rest-frame wavelength A. Although the extinction may make 
sources harder to detect, the dependence on wavelength is ex- 
pected to be smooth enough so as not to distort the signature 
of H I absorption. Specifically, near the Lya wavelength, a 
given neutral column produces extinction 



A(A) w 1.7 



1-0.013 



X— X a 
10 A 



10 21 cm" 



(16) 



where all wavelengths in this expression are measured in the 
rest frame. Note that this expression assumes the Galactic 
dust-to-gas ratio, while a high-redshift galaxy is expected to 
have a lower mean metallicity and dust content (although high 
metallicity may be found in star-forming regions due to strong 
local metal enrichment). 

3. THE BASIC MODEL: THREE EXAMPLES 

In this subsection we illustrate our predictions for several 
different realistic sets of parameters. We consider a GRB 
and a quasar, each at redshift 7 and having formed in a re- 
gion that had not yet been reionized by other sources (or had 
recombined after a possible earlier reionization phase; see 
§1). The quasar is assumed to have an observed flux F\ of 
10" 17 erg cm" 2 s" 1 A (typical of the SDSS high-redshift 
quasars) and to have been radiating for ts = 2 x 10 7 yr a t the 
Eddington luminosity, which implies in our model (see jj2.3> 
a 1 .7 x 10 9 M Q black hole at the center of a halo of total mass 
4.7 x 10 12 M Q . In the other example, we assume that the GRB 
goes off in a halo of mass M = 4 x 10 8 M Q (typical of z = 7 
GRBs; see ^3}, which hosts a galaxy of stars; we consider ei- 
ther a normal IMF with / esc = 10% and t s = 1 x 10 7 yr, or a 
Pop III IMF with / esc = 90% and t s = 1 x 10 s yr. 



Figure ^ shows the resulting absorption profiles in each of 
these cases, and also presents the physical properties of the 
gas within the H II region. A useful parameter is the apparent 
line-of-sight position AD relative to the source redshift zs, 
derived directly from the observed wavelength by assuming 
Hubble flow. The conversion, for a small wavelength differ- 
ence AA relative to Lya at zs, is 



AD = 0.370 



AA \ 

Tok) 



1+zs 



-5/2 



0.136 



-1/2 



Mpc 



(17) 

The figure shows that the contributions to r from the damping 
wing (rdamp) and from resonant absorption (t ks ) are compara- 
ble. Both become very large as the edge of the H II region is 
approached on the blue side (at AD < 0). On the red side, 
however, the smooth decline of Td amp with AD contrasts with 
the sharp drop in r res . Clearly, regardless of the presence or 
ab sence of t he damping wing, the signature of infall identified 
by Barkana & Loeb ( 2003) is expected, i.e., a sudden and sig- 
nificant drop in the transmitted flux at the wavelength where 
resonant absorption starts to operate. Note that in order to cal- 
culate the damping wing we assume a neutral IGM along the 
line of sight down to z = 6.5, although Td amp is very insensitive 
to the precise lower limit [see eq. Q]. 




-4 -2 o 
AD [Mpc] 



-0.2 -0.05 0.05 
AD [Mpc] AD [Mpc] 



FIG. 1 . — Properties of the gas in the H II region versus relative apparent 
position AD. Panels correspond to a GRB host with a normal stellar IMF 
(right panels) or a Pop III IMF (middle panels), and to a quasar (left panels), 
all at redshift 7 in a region that had not been pre-ionized (note the different 
AD ranges; see text for the values of the input parameters). From bottom to 
top, panels show the Lya optical depth due to resonant absorption (t ics ; solid 
curves) and due to the damping wing (r,j am p; dashed curves); the clumping 
factors CLya (solid curves) and C VHI (dashed curves); the gas overdensity 
due to infall; and the actual distance of each gas element from the ionizing 
source (solid curves), as well as the distance that the same gas element would 
have reached in the presence of pure Hubble expansion with no infall (dashed 
curves). 

The figure shows the strong effect of infall, which produces 
resonant absorption well into the red wing (i.e., at AD > 0). 
In the case of pure Hubble flow, r res would drop to zero right 
at AD = 0, and the drop would be gradual because of the in- 
creasing ionizing intensity at small distances from the source. 
When infall is included as part of a more realistic model, how- 
ever, the infall velocity results in absorption of photons at 
Az > 0; furthermore, the value of r res does not show a grad- 
ual drop since the high infall overdensity (and, therefore, high 
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recombination rate) compensates for the high ionizing inten- 
sity felt by gas near the source. The sudden drop in r les is 
due to the inclusion of an accretion shock, which fully ion- 
izes the gas and cuts off absorption very close to the source. 
The inclusion of a clumping distribution, and the dominance 
of voids in allowing transmitted flux, leads to a strong reduc- 
tion in r res (i.e., C^ya is significantly smaller than C XH1 ). In the 
quasar and Pop III star cases, helium is fully ionized through- 
out the H II region, but in the normal IMF case, the GRB host 
galaxy produces only a small He III region; the optical depth 
jumps at this radius, due to the lower gas temperature in the 
He II region, and the transition can be seen at AD = 0.02 Mpc. 

The effective clumping factor for transmission (CL yQ ) varies 
over the range 0.03-1 . On the other hand, the clumping factor 
C Xm is always much larger, although its value (~0.5-3) is still 
significantly lower than the values > 10 that previous studies 
have used which did not account for the rarity of high clump 
densities along a given line of sight. Since it is defined in 
terms of a ratio of neutral fractions [see eq. 0], C Xm is essen- 
tially independent of both A and the ionizing intensity, and is 
simply a function of A max . For the cases shown in the figure, 
the proper radius of the H II region and the value of A max are 
respectively 0.063 Mpc and 1.8 (right panels), 0.28 Mpc and 
3.7 (middle panels) and 4.3 Mpc and 31 (left panels). On the 
other hand, CL ya is dominated by low values of Ac, is nearly 
independent of A max , and is essentially a function of the no- 
clumping (C = 1) optical depth, which itself is proportional to 
A 2 R , where we used the fact that the source ionizing inten- 
sity (which affects the neutral hydrogen fraction) is inversely 
proportional to R 2 . Indeed, when the C = 1 optical depth is 
large, the transmission in the voids (where r res oc A c is low) 
becomes increasingly dominant in the calculation of the over- 
all mean transmission. Note that our results dif fer quanita- 
tively from the clumping model of Haiman ( 2002), who found 
an effective clumping factor for transmission of ~ 0.04 when 
T res = 1, while we obtain CL yQ ~ 0.4 in this case. 

The figure also shows that A ~ 20 at the shock radius, al- 
though the infall pattern depends on the halo mass and red- 
shift; at a given redshift, higher-mass halos correspond to 
rarer density peaks, which tends to strengthen the infall pat- 
tern around them, but in standard models there is less power 
on large scales and this works in the opposite direction [see 
Barkana (2004)]. Gas with large A is located at a real radius 
R that is much smaller than the radius the same gas element 
would have reached in the absence of infall; this is consistent 
with mass conservation. Observations of r res can in princi- 
ple probe the density profile of cosmological gas infall all the 
way from the shock radius out to the radius of the H II region, 
a range that represents ~ 1 .5 orders of magnitude of R. Note 
that for each gas element, the total velocity relative to zs is 
H(z s )lS.D, while the Hubble velocity component is H(zs)R, 
where the Hubble constant at high redshift is accurately ap- 
proximated by 

The main conclusion from this figure is that a quasar af- 
fects a much larger surrounding region than a typical GRB 
host, both in terms of infall and ionization. Specifically, the 
region of significant infall (measured by A > 2) is larger by 
a factor of ~ 10 in radius, or 1000 in volume. Likewise, the 
H II region radius is larger by a factor > 10 even compared 
to a GRB host with an extreme Pop III population of stars. 
Furthermore, within the H II region itself, r res is of order unity 



for a quasar, and of order 10 in the GRB case. This means 
that absorption in this region can be easily probed in quasars, 
while for GRB hosts the absorption is near total except well 
into the damping wing; note that at larger AD, beyond the 
limits shown in the figure, Td amp continues to decrease roughly 
inversely with AD- ADn, where ADn < is the apparent po- 
sition of the edge of the H II region. Thus, cosmological infall 
is pronounced around quasars and its effects are observable 
over a large region, while it is weak in GRB hosts and essen- 
tially unobservable. On the other hand, the damping wing can 
be observed in GRB hosts over a large region to the red side 
of resonant absorption, while in quasars Td amp <C 1 except in 
the region where r res is of comparable importance. 

Note that we have not solved for exact radiative transfer 
in the expanding H II region, but instead have assumed that 
the region is optically thin when we determined the ioniz- 
ing intensity at various positions within the region. In or- 
der to check this assumption, we calculate the optical depth 
Ti on to hydrogen-ionizing photons from the source being ab- 
sorbed as they travel out to the edge of the H II region. This 
optical depth depends on the intensity and gas density at all 
radii through the region. In order to place an upper limit on 
the optical depth, we evaluate it for photons with an energy 
just above the ionization threshold of hydrogen. For the three 
cases shown in the figure, T; on thus defined equals 0.019 (nor- 
mal GRB host), 0.027 (Pop III GRB host) and 0.015 (quasar). 
The fact that T; on <C 1 justifies the neglect of radiative absorp- 
tion and scattering in our calculations. 

4. DETECTING COSMOLOGICAL INFALL AROUND QUASARS 

As illustrated in the previous section, observations of Lya 
absorption in quasars can probe the cosmological infall pat- 
tern induced by the dark matter halo surrounding each quasar. 
In this section, we predict the dependence of the absorption 
pattern on the quasar redshift and flux. If these trends pre- 
dicted by theory can be checked directly in observations, then 
this would test the fundamental properties of quasars and of 
dark matter halos, as well as the correlation between them. 

Quasars have broad-band spectra that, due to their high lu- 
minosity, can be measured accurately even at high redshift. 
However, unlike GRB afterglows, the spectra are composed of 
emission lines broadened by multiple components of moving 
gas, and the difficulty in predicting the precise shape of, e.g., 
the Lya line in any particular quasar presents a significant sys- 
tematic error in determinations of the absorption profile. A 
promising approach is to adopt the parametrized form of the 
emission profile that best fits the line shape of most quasars 
at low redshift, and then fix the parameters for each quasar 
based on th e part of the Lya line th at does not suffer resonant 
absorption ( Bar kana & LoebT20 03). Quasars also have impor- 
tant advantages. The great intensity of quasars highly ionizes 
their host galaxies and reduces the amounts of gas and dust, 
thus simplifying the analysis. Furthermore, if the relatively 
high metallicity found for high-redshift quasars is typical of 
quasars at even higher redshifts, then various metal lines may 
be detected in addition to Lya, allowing an independent mea- 
sure of the host redshift and of the velocity structure of the 
gas producing the broad emission lines. 

We can easily understand the rough scaling with the ob- 
servable parameters (i.e., quasar luminosity and redshift) of 
the properties of the red absorption drop due to infall. The 
luminosity of a quasar at redshift zs is proportional to Mbh 
times /g, the fraction of the Eddington luminosity that is ac- 
tually emitted by a given quasar. We assume here a general 
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power-law relation Mbh oc V" between the black hole mass 
and the halo circular velocity (note that we adopt n = 5 in the 
numerical results, as discussed in 32. 3> . The distance to the 
accretion shock is proportional to the halo virial radius, and 
the infall velocity of the pre-shock gas is proportional to the 
halo circular velocity. Thus, the wavelength offset of the ab- 
sorption drop (measured redward from the central wavelength 
of the source Lya emission) is 

/ L \ 1/n 

AA drop cx(l+z s ) ( — J . (19) 

The halo mass is given by M oc (1 +zs)~ 3 ^ 2 V e 3 , or 

/ L \ 3/ " 

Moc(l+z s )~ y2 I— J . (20) 

The optical depth at the absorption drop is the Gunn-Peterson 
optical depth at that redshift, times the relative gas density 
A, times the neutral hydrogen fraction. The neutral fraction 
is proportional to the gas density over the ionizing intensity, 
where the intensity is proportional to the quasar luminosity 
over the square of the distance to the accretion radius. Thus, 
the optical depth scales as 

T (2-n)/n 

r drop ^A 2 (l +zs) y2 ^j^- , (21) 
Je 

where A depends on M and z through the infall model. Note 
that the dependence of optical depth directly on halo mass is 

r drop oc A 2 ( 1 + zsf-" >/2 fE W 2 ^ 3 . (22) 

Also note that clumping f § 12.21 modifies slightly these scal- 
ings of the optical depth. Observing the absorption drops 
in quasars over a wide range of luminosity and redshift will 
probe the population of quasars and the properties of their ha- 
los through the values of fg, A, and n. 

Figure[2]shows the fraction of the quasar's flux that is trans- 
mitted through the surrounding, infalling IGM. We assume a 
quasar at redshift 6 that had turned on after reionization, and 
we therefore include only resonant absorption. We consider 
a quasar with a flux typical of the recent SDSS high-redshift 
quasars, and also compare to quasars with fluxes lower by a 
factor of 10 or 100, applicable to future surveys. In each case, 
we translate the flux to black hole mass by assuming that the 
quasar radiates at the Eddington luminosity. We also compare 
to the case where the quasar is instead assumed to be radiating 
at only 10% of Eddington; in this case the black hole must be 
more massive in order to produce a given observed flux, im- 
plying a larger surrounding halo, and an absorption drop more 
strongly shifted to the red (corresponding to a larger infall ve- 
locity). As implied by the scalings presented above, fainter 
quasars at a given redshift have a smaller incursion of absorp- 
tion into the red wing of Lya; they also have a larger r^op 
since the slightly smaller A around lower-mass hosts does not 
compensate for the lower luminosity of the quasars. 

The transmitted flux fraction is also shown in Figure[3] but 
for various quasar redshifts. For a given observed flux, lower- 
redshift quasars are situated inside smaller hosts, and the ab- 
sorption profiles are seen to have a smaller AAd r0 p along with 
a slightly lower T± op . For zs = 7, a quasar emitting inside a 
still neutral region is distinguished by the appearance of the 
damping wing. In addition, such a quasar has a significantly 
smaller red absorption drop because of the higher tempera- 
ture of gas reionized by the relatively hard quasar spectrum; 
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FIG. 2. — Transmission fraction of a quasar's flux through the infalling 
intergalactic medium, as a function of wavelength (or velocity) relative to 
the center of the quasar's emission line. The quasar is assumed to have an 
observed flux F\ of 1 (solid curves), 0.1 (dashed curves), or 0.01 (dotted 
curves), in units of 10~ 17 erg cm -2 s A . In each pair of curves, the 
quasar is assumed to be radiating at its Eddington luminosity (top) or at 10% 
of Eddington (bottom). All curves are at redshift 6 and assume that the quasar 
had turned on after reionization. 
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FIG. 3. — Transmission fraction of a quasar's flux through the infalling 
intergalactic medium, as a function of wavelength relative to the center of the 
quasar's emission line. The quasar is assumed to have an observed flux F\ of 
10~ 17 erg cnT 2 s _1 A -I , and to be in a post-reionization universe at redshift 
7 (solid curves), 5 (short-dashed curves), or 3 (dotted curves). Also shown is 
the same redshift 7 quasar but assuming that it reionized its surroundings for 
the first time over a time ts = 2 X 10 7 yr (long-dashed curves). In each pair 
of curves, the quasar is assumed to be radiating at its Eddington luminosity 
(top) or at 10% of Eddington (bottom). 



the recombination rate is lower in such gas, thus reducing the 
hydrogen neutral fraction and the Lya optical depth. 

In both figures, apparent at large distances from the quasar 
is the contribution of the mean ionizing background, which is 
included at redshift 6 and below but is most prominent at red- 
shift 3. In addition, the appearance of the red-most absorption 
feature, and the shape of the rise and then the fall of the trans- 
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mission toward shorter wavelengths, vary widely depending 
on the properties of the quasar, its halo, and the surrounding 
infall region. In principle, the precise shape of the absorp- 
tion profile can be used to measure the detailed gas density 
profile in addition to the basic parameters of the quasar and 
its host halo. However, these transmission profiles cannot be 
directly observed; for any particular quasar, the transmission 
profile operates on the intrinsic profile of the quasar emission, 
and density fluctuations along the line of sight change the 
transmission relative to the ensemble-averaged profile. While 
these complications can be modeled, we have focused here 
on the observable features that are independent of these un- 
certainties and that are recognizable even in individual spec- 
tra, although the precise values may fluctuate compared to 
the ensemble-averaged ones; specifically, these features are 
the magnitude of the red absorption drop, and its location in 
wavelength, which is most useful if the quasar redshift is mea- 
sured accurately, e.g., from other emission lines. For a sta- 
tistical sample of quasars, our results should best match the 
ensemble-averaged spectrum. 

5. PROBING REIONIZATION WITH 7-RAY BURST AFTERGLOWS 

Gamma-Ray Burst (GRB) explosions from the first gen- 
eration of stars offer a particularly promising opportunity 
to probe the epoch of reionization. Young (days to weeks 
old) GRBs are known to outshine their host galaxies in 
the optical regime at z ~ 1-3. In hierarchical galaxy for- 
mation in ACDM, the characteristic mass and hence op- 
tical luminosity of galaxies and quasars declines with in- 
creasing redshift; hence, GRBs should become easier to ob- 
serve than galaxies or quasars at in creasing redshift. If the 
first s t ars are pre ferentiall y massi ve ( Abel, Brvan, & Norman 
120021 iBrom m. Copp i. & Larsonl 120021) and GRBs originate 
during the formation of compact stellar remnants such as neu- 
tron stars and black holes, then the fraction of all stars that end 
up as GRB progenitors may be even higher at early cosmic 
times. 

Similarly to quasars, GRB afterglows possess broad-band 
spectra which extend into the rest-frame UV and can probe 
the ionization state and metallicit y of the IGM out to the 
epoch when it was reionized (Lamb & Reichart 2000). Sim- 
ple scaling of the long-wavelength spectra and temporal evo- 
lution of afterglows with redshift implies that at a fixed time 
lag after the GRB trigger in the observer's frame, there is 
only a mild change in the observed flux at i nfrared or radio 
wavel engths as the GRB redshift increases ( Ciardi & Loeb 
2000). Hence, GRBs provide exceptional lighthouses for 
probing the universe at high redshift. The infrared spectrum 
of GRB afterglows occurring around the reionization redshift 
could be taken with future telescopes such as the James Webb 
Space Telescope (JWST; see http://www.jwst.nasa.gov/!, as 
a follow-up on their early X-ray localization with the Swift 
satellite. In this section, we explore how such a spectrum 
could reveal the damping wing due to neutral IGM, and thus 
probe reionization. 

If GRBs originate from stellar remnants of high-mass stars, 
then they are expected to occur in host galaxies with a proba- 
bility proportional to the SFR of each galaxy. In order to spec- 
ify a typical host, we estimate as follows the SFR of galax- 
ies in halos of various mass M at redshift z = 7. We assume 
that the SFR in a given halo is proportional to M divided by 
a star-formation timescale which is determined by the halo 
merger history. Specifically, we estimate the age of gas in a 
given halo using the average rate of mergers which built up 



the halo. Based on t he extended Press-Schechter formalism 
JLacev & Coldl 19931) . for a halo of mass M at redshift z, the 
fraction of the halo mass which by some higher redshift zi had 
already accumulated in halos with galaxies is 



t? / \ f f 1-69/ D(Z2)- 1-69/ D(z) \ 



(23) 



where D(z) is the linear growth factor at redshift z, S(M) is 
the variance on mass scale M (defined using the linearly- 
extrapolated power spectrum at z = 0), and M m ; n (z2) is the 
minimum halo mass for hosting a galaxy at zi- We assume 
that prior to the final phase of reionization this minimum mass 
is determined by the minimum virial temperature of ~ 10 4 K 
required for efficient atomic cooling in gas of primordial com- 
position. We estimate the typical age of gas in the halo as the 
time since redshift zi where Fm(z,Z2) = 0.5, so that 50% of the 
gas in the halo has fallen into galaxies only since then. The 
result at z = 7 is that M min (z = 7) = 1 .3 x 10 8 M Q , 50% of the 
total SFR occurs in halos up to M = 4 x 10 8 M Q , and 90% in 
halos up to M = 5 x 10 9 M Q . 

The GRB afterglow spectrum can be modeled as a set 
of power-law segments. Regardless of the particular power 
law which applies to the Lya region, since the wavelength 
band relevant for H I absorption is narrow, the intrinsic after- 
glow emission is essentially constant over this region; e.g., 
at z = 7, the flux level should vary by < 1 % over a range 
of A A = 100 A around the redshifted Lya wavelength. Thus, 
the GRB afterglow offers a unique opportunity to detect the 
signature of IGM absorption, if the afterglow spectrum is 
at all similar to the theoretically expected smooth spectrum. 
The followin g numerical expressions and figu res are based on 
the model of Sari, Piran, and Naravan ( 1998), corrected for a 
jet geometry us ing Sari, Piran, and Halpern ( 1999) [see also 
Waxman & Gruzinovf ll 9991) 1. 

Based on modeling of observed afterglows, we assume typ- 
ical afterglow parameter s of E = 1 51 erg for the total (not 
isotropic) energy ( Frail et al, 2001), a jet opening angle of 
9 = 0.07 rad, an external shock occurring in a surrounding 
medium of uniform density n = 1 cm" 3 , electrons acceler- 
ated to a power-law distribution of Lorentz factors with index 
p = 2.4, and a fraction e B = 10% of the shock energy going to 
the magnetic field and e e = 10% to the accelerated electrons. 
We assume adiabatic evolution, and neglect synchrotron self- 
absorption which is only important at much lower frequencies 
than we consider. 

Figure|4f shows quantitatively the detectability of GRB af- 
terglows at redshifts z = 5-15. The figure indicates that the 
redshifted Lya region of GRB afterglows remains sufficiently 
bright for a high-resolution study with JWST up to several 
weeks after the GRB explosion, throughout the redshift range 
we consider. Each light-curve of afterglow flux shows a cou- 
ple of breaks over the range of times considered in the figure. 
We let v m and v c be the characteristic synchrotron frequencies 
(redshifted to the observer) of, respectively, the lowest-energy 
electron in the accelerated power-law distribution, and the 
lowest-energy electron that radiates efficiently enough to lose 
most of its energy. At early times the afterglow is in the fast 
cooling regime (y m > u c ), when all the accelerated electrons 
cool efficiently, and the observed frequency v a (correspond- 
ing to redshifted Lya) is intermediate between the two char- 
acteristic frequencies, leading to a spectral shape F v oc 5 

2 Detection thresholds in Figure|4]were obtained using the JWST calcula- 
tor at http://www.stsci.edu/jwst/science/jms/index.html . 
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near Lya. The first break occurs when v a rises above v m at 



fobs - 6.8 



1/3 



10 51 erg/ V°- 07rad / Vo1 



-2/3 



1/3 



(24) 



After this break, ^ Q is bigger than both !/ m and i/ c , and F v oc 
i/" p / 2 near Lya. The second break, which is due to sideways 
expansion of the jet, affects the light-curve but not the spectral 
slope near Lya. It occurs at 



fobs =3.5 



1/3 



10 51 erg/ V°- 07rad 

n \-i/3 fl+z s . , 

days 



1 cm 



(25) 




FIG. 4. — Detectability of high-redshift GRB afterglows as a function of 
time since the GRB explosion as measured by the observer. The GRB af- 
terglow flux is shown (solid curves) at the redshifted Lya wavelength. Also 
shown (dotted curves) is the detection threshold for JWST assuming a spec- 
tral resolution R = 5000 with the near infrared spectrometer, a signal to noise 
ratio of 5 per spectral resolution element, and an exposure time equal to 20% 
of the time since the GRB explosion. In each set of curves, a sequence of 
redshifts is used, z = 5, 7, 9, 11, 13, and 15, respectively, from top to bottom. 

As noted earlier, afterglows fade only slowly with redshift, 
since a given observed time corresponds to an earlier source 
time for high redshift afterglows, which are thus intrinsically 
brighter. Indeed, after the second break, the Lya flux at a 
given f b s is enhanced by this effect in proportion to (1 +zs) 
to the power p, a factor that for p = 2.4 almost overcomes the 
fading according to the inverse squared luminosity distance. 
GRB afterglows may be even more easily detectable than in- 
dicated in figure 0] if a significant fraction exhibit an early 
UV-optical flash from a reverse external shock; e.g., the flash 
in GRB 990123 at z = 1.6 was detected at a flux level of — 1 
Jy at r est-frame A ~ 2000A about 50 seconds after the burst 
began jAkerlofet all 19991) . 

Relative to other possible sources, GRB afterglows offer 
the advantage of having smooth, power-law spectra, which 
helps GRBs avoid the systematic uncertainties associated with 
the need to model the Lya emission profiles of galaxies and 
quasars. In addition, GRBs are expected to occur mainly in 



relatively low mass halos, which induce only weak infall of 
the surrounding gas and produce insignificant H II regions, 
as illustrated in § [3] Metal lines from the host may be seen 
in absorption in the afterglow spectrum, thus allowing an ac- 
curate redshift measurement. Emission from the host galaxy 
should not significantly contaminate the early afterglow, but it 
may make the galaxy marginally detectable with JWST once 
the afterglow has faded away. For instance, 10% of the SFR 
at z = 7 occurs in halos above M = 5 x 1O 9 M , and a halo 
of this mass is expected to have a SFR of ~ 1M© per year, 
which translates [for a normal IMF] to a 1500 A continuum 
flux density of ~ 15 nJy, and a flux density in the ~ 10A band- 
width of the Lya line emission from the galaxy of ~ 500 nJy. 
If the star-forming regions where GRBs are expected to oc- 
cur contain significant column densities of gas and dust, then 
they may c ontam inate the signal from IGM absorption, as dis- 
cussed in § 12.61 

Figure [5] illustrates GRB afterglow absorption profiles (top 
panel) and explores the possibility of contamination from dust 
and neutral gas in the host (bottom panel). If absorption from 
the host is negligible then GRB afterglows offer a clear and 
direct measurement of absorption by the damping wing of H I 
in the IGM. In particular, the sharp cutoff expected if only res- 
onant absorption is included (bottom panel, dot-long dashed 
curve) is easily distinguishable from the gradual cutoff char- 
acteristic of absorption by the damping wing (bottom panel, 
solid curve). In this case, the IMF can be constrained as well, 
since a strongly-ionizing IMF produces a gradual cutoff as 
well as a sharp drop due to resonant absorption (top panel, 
dashed curves); this occurs since the neutral gas lies further 
away (in the blue wing), producing a relatively weak damp- 
ing wing that is interrupted by the ubiquitous red-wing drop 
due to resonant absorption. 

However, as illustrated in the bottom panel of figure[5] con- 
tamination by absorption in the host galaxy may hamper stud- 
ies of reionization and of the IMF based on absorption pro- 
files. The three lower curves include extinction as might be 
produced by an H I column density of 1 2 1 cm" 2 [see eq. J16i 1 . 
Because of its relatively smooth spectrum, extinction does 
not significantly change the shape of the absorption profile; 
moreover, it is possible to take advantage of the broad-band 
afterglow spectrum in order to derive the effect of extinction 
near redshifted Lya by measuring the extinction over a much 
wider range of wavelengths than shown in the figure and inter- 
polating into the critical bandpass. However, extinction does 
lower the overall flux and thus increases the integration time 
required to measure an accurate spectrum. 

As shown in figure |4] at f ODS = 1 day, a minimum flux of 
F v ~ 10~ 6 Jy is needed for a high-resolution spectrum with 
JWST. Assuming a spectrum obtained with this sensitivity, 
we consider whether absorption by a large H I column den- 
sity in the host could be confused with the signature of the 
IGM damping wing before reionization. We consider only 
the F v > 10~ 6 Jy portion of the profile with extinction and 
full absorption due to the IGM (figure|5j bottom panel, dotted 
curve), and attempt to reproduce this profile by replacing the 
IGM damping wing with a damped Lya system in the galac- 
tic host. If the host redshift is precisely known, the different 
shapes of the two gradual cutoffs preclude any possibility of 
confusion between them. For example, a damped Lya system 
with /Vhi = 2 x 10 2l) cm" 2 (short-dashed curve) is 1.5 orders 
of magnitude too low at AA ~ 10A and 10-20% too high at 
AA = 30-100A. 
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FIG. 5. — IGM absorption profiles in GRB afterglows, presented in terms of 
the flux density F u versus relative observed wavelength AA. All curves adopt 
z = 7 (assumed prior to the final reionization phase) and the typical halo mass 
M = 4 X 10 8 Mq expected for GRB host halos. Top panel: Predicted spectrum 
including IGM H I absorption (both resonant and damping wing), for host 
galaxies with t$ = 10 7 yr, / esc = 10% and a normal IMF (solid curves) or t$ = 
10 s yr, /esc = 90% and a Pop III IMF (dashed curves). The observed time after 
the burst is one hour, one day, and ten days, from top to bottom, respectively. 
Bottom panel: Predicted spectra one day after a GRB for a host galaxy with 
t$ = 10 7 yr, /esc = 10% and a normal IMF. Shown is the unabsorbed GRB 
afterglow (dot-short dashed curve, essentially horizontal), the afterglow with 
resonant IGM absorption only (dot-long dashed curve), and the afterglow 
with full (resonant and damping wing) IGM absorption (solid curve). Also 
shown, with 1.7 magnitudes of extinction, are the afterglow with full IGM 
absorption (dotted curve), and attempts to fit this profile with a damped Lya 
absorption system in the host galaxy (dashed curves; see text). 



However, if the source redshift is not measured exactly then 
a much closer match is possible. The damped Lya system 
(long-dashed curve) that most closely simulates the effect of 
the IGM damping wing in this case has Nm = 7.7 x 10 20 cm" 2 
and is shifted blueward by AA = -15.4A, equivalent to Az 
= 0.013 or a velocity blueshift of 500 km s" 1 . This match, 
although it appears fairly close in the figure, is off at vari- 
ous wavelengths by a difference ranging between -6% and 
+6%, and this difference should be easily measurable. How- 
ever, the relatively close match in this comparison of a pure 
IGM damping wing with a pure damped Lya system suggest 
that in a realistic situation, when both may contribute, this 
partial degeneracy will limit the constraints that can be de- 
rived on the various parameters of interest. Note, though, 
that since the halo circular velocity is only ~ 24 km s -1 , 
the peculiar velocity of the GRB (which may arise from the 
circular velocity of the disk) is negligible, and an accurate 
redshift measurement can help to bre ak the near-degeneracy . 
Regarding another possible difficulty, Miralda-Es cudel (119981) 
showed that ~ 10% uncertainties in the cosmological pa- 
rameters would contribute further to the near-degeneracy be- 
tween a damped Lya system and the IGM damping wing, but 
this factor continually becomes less important as various cos- 
molo gical probes yield in creasingly precise parameter values 
(e.g JSpergel et all2003l) . 

The hydrogen absorption within the host galaxy of the 
GRB should be reduced b y the early ionizing r adiation of 
the GRB afterglow itself JPerna & Loebl fl998). The re- 



combination time of the surrounding interstellar medium (~ 
lO 5 ^ 1 yr) is long compared to the duration of the after- 
glow emission, where no is the medium density in units of 
cm -3 . Hence, each ionizing photon will eliminate a neu- 
tral hydrogen atom from the interstellar medium along the 
line of sight. Given that the number of hydrogen atoms 
out to a distance r is oc nor 3 , the radius ionized by a given 

afterglow varies as n . iDraine & Had J2002I) showed 
that an early optical/UV flash should ionize hydrogen and 
destroy dust and molecular hydrogen out to a distance of 
~ 10 19 cm in surroundings of no ~ 10 3 — 10 4 (see also 
i Perna & Loeblll998l IWaxman & Draine! 120001: lOraind [ 20001: 
Fruchter. Krolik. & Rhoadsl 1200 It iPerna & Lazzatil booi 
Perna, Lazzati, & Fiore 2003). Therefore, the early afterglow 
emission could significantly reduce the damped Lya absorp- 
tion within a host molecular cloud {no ~ 10 4 , r ~ 10 pc) or a 
line of sight passing at an angle through a host galactic disk 
(no ~ 1, r ~ 100 pc). However, the afterglow's characteristic 
sphere of influence is negligible compared to the intergalactic 
scales of interest for reionizing the IGM. 

6. CONCLUSIONS 

We have calculated the profile of Lya absorption for GRBs 
and quasars during and after the epoch of reionization, includ- 
ing the effects of ionization by the source and its host galaxy, 
clumping, and cosmological infall of the surrounding cosmic 
gas. We find that GRBs are the optimal probes of the ion- 
ization state of the surrounding IGM as they are most likely 
to reside in dwarf galaxies that do not significantly perturb 
their cosmic environment. On the other hand, bright quasars 
reside in the most massive galaxies and probe the strong grav- 
itational effects of these galaxies and their halos on the IGM. 

In order to isolate the intergalactic contribution to the Lya 
damping wing, it is necessary to measure the source redshift 
accurately. Here again, GRB afterglows are simpler to in- 
terpret. In quasars, although there are often strong emission 
lines in addition to Lya, bulk motions can shift them substan- 
tially from the systemic redshift of the host galaxy. For exam- 
ple, some broad emission lines such as C IV are typically off- 
set by ~ 1000 km s , but low-ionizatio n broad lines such a s 
Mg II are only offset by ~ 200 km s" 1 jRichards et alJl2002l) . 
The Mg II line can be detected at high redshift in the near- 
infrared and, indeed, it has been observed for the high est red- 
shift quasar known JWillott. McLure. & Jarvisl 120031) . Red- 
shifts accurate to better than 100 km s" 1 may be obtainable by 
detecting narrow lines such as [O III], perhaps marginally pos- 
sible from the ground or else doable with JWST. For either a 
GRB or a quasar, however, perhaps the best method for deter- 
mining its redshift is to detect molecular emission or absorp- 
tion lines in the interstellar medium of its host galaxy. The re- 
sulting measurement is likely to be accurate to within the disk 
velocity dispersion, which should be rather low, particularly 
in the low-mass hosts that are typical of GRBs. These lines 
can be detected at high redshift, as demonstrated by the CO 
emissio n from the host ga laxy of the highest-redshift quasar 
known ( Fabian et al. 2003). Note that another possible dis- 
advantage of quasars is that their massive hosts may develop 
significant peculiar velocities as they are likely to reside in 
rare high-density environments such as groups of galaxies. 

Supernova and quasar outflows could change our predic- 
tions if they reached out from the galaxy far enough to affect 
the gas at the distances where it produces the absorption in 
our model. Although outflowing gas can be detected spec- 
trally and the outflow velocity deduced, it is generally very 
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difficult to determine the spati al scale of an o utflow around a 
high-redshift galaxy. Adelbergeret al. (2003) have detected 
a statistical signature based on six objects that suggests that 
Lyman break galaxies at z ~ 3-4- produce outflows that drive 
away surrounding gas out to a radius of ~ 0.5 h~ l comov- 
ing Mpc. If this preliminary signal is correct, and if these 
Lyman break galaxies lie in M > 10 12 M Q host halos, then 
these outflows reach a distance comparable to the halo ac- 
cretion shock and may begin to marginally affect our model 
predictions. Note that direct absorption by the outflowing gas 
should generally be distinguishable from the absorption pat- 
terns that we predict. Absorption from outflows should occur 
in the blue wing and thus not interfere with the red absorption 
feature due to infall. Furthermore, outflowing gas is likely to 
be highly enriched and can be probed (or ruled out) in each 
object using metal lines. 

High-redshift GRB afterglows are likely to open a new 
frontier in extragalactic astronomy ("GRB cosmology'), by 
providing important new clues about their progenitors — the 
first stars — and about the reionization history of the interven- 



ing IGM. Afterglow spectra can also be used to explore the 
metal enrichment history of the IGM thr ough the detection 
of in tergalactic metal absorption lines ( Furlan etto & Loebl 
120031) . The comparison between observed spectra of high- 
redshift afterglows and detailed cosmological simulations of 
the IGM (involving hydrodynamics and radiative transfer dur- 
ing the epoch of reionization) will provide a better under- 
standing of the process of galaxy formation in the young uni- 
verse, only ~ 10 s — 10 9 years after the big bang. 
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